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ABSTRACT 

The frequency of short gamma-ray bursts (GRBs) in galaxies with distinct star formation histories 
can be used to constrain the lifetime of the progenitor systems. As an illustration, we consider here 
the constraints that can be derived from separating the host galaxies into early and late types. On 
average, early-type galaxies have their stars formed earlier than late-type galaxies, and this difference, 
together with the time delay between progenitor formation and short GRB outburst, leads to different 
burst rates in the two types of hosts. Presently available data suggest, but not yet prove, that the 
local short GRB rate in early-type galaxies may be comparable to that in late-type galaxies. This 
suggests that, unlike Type la supernovae, at least half of the short GRB progenitors that can outburst 
within a Hubble time have lifetimes greater than about 7 Gyr. Models of the probability distribution 
of time delays, here parametrized as P(t) oc t™, with n > — 1 are favored. This apparent long time 
delay and the fact that early-type galaxies in clusters make a substantial contribution to the local 
stellar mass inventory can explain the observed preponderance of short GRBs in galaxy clusters. 

Subject headings: gamma rays: bursts - stars:formation - cosmology: observations - galaxies:formation 



1. INTRODUCTION 

The progenitors of short duration, hard spectrum, 
gamma-ray bursts (GRBs) are not yet well identified. 
Even with the recent localizations of a handful of short- 
hard GRBs flB i oorri^t~aTl 120061: iGefirels et all 120051: 
Villasenor et al.l 120051: IBerger et all 120051 : I Fox et all 
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no transient 

emission has been found that directly constrains the na- 
ture of the progenitor system. The current view of most 
researchers is that GRBs arise in a very small fraction of 
stars that undergo a catastrophic ene rgy release event to- 
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is dedicated to understanding the different progenitor 
scenarios and trying to determine how the progenitor 
and the burst environment can affect the observable 
burst and afterglow characteristics (e.g., iLee et all 
2005). The lifetime of progenitors of short bursts can 
be meaningful ly constrained by properties of their host 
galaxies (e.g., iFox et aT1l2lM iGaT-Yam et alJl2005h . We 
suggest here a method to constrain the lifetime of the 
progenitors by making use of the star formation histories 
of their host galaxies. 

Based on th e short burst afterglows localized so far 
(|Berger I [20 07). two of the three relatively nearby (z < 
0.3) events (GRB 050509b and GRB 050724) are plau- 
sibly associated w ith galaxies exhib i ting characteristic 
early- type spectra (|Berger et al.ll2005l : lBloom et al.ll2006t 

1 Institute for Advanced Study, Einstein Drive, Princeton, NJ 
08540; zhengz, enrico@ias.edu 

2 Hubble Fellow 

3 Chandra Fellow 



iGehrels et al]|2005b [Prochaska et al.ll2006f h In the other 
case (GRB 050709), the host galaxy shows signature 
of a dominant stel lar population with age of ~ lGyr 
(|Covino et al.ll2006f) and it also exhibits st rong emission 
lines t hat indica t e ong o ing star formation (ICovino et all 
120061: IFox et alJ 120051: iProchaska et all [2006J. There 
is also independent support that, at least t wo o f the 
nine short bursts localized so far (Berger 20Q 3) are 
associated with clus t ers of galaxies ( Bloom et al.ll2006l 
IGladders et~aT]|2005t iPedersen et al.ll2005l ). In contrast 



to what is found for long-soft GRBs, for which all of the 
confirmed host galax i es are actively forming stars (e.g., 
iTrentham et~aLll2002l : IChristensen et al.ll2004h . these ob- 
servations clearly signify that, like Type la supernovae, 
short GRBs are triggered in galaxies of all types. What is 
more, it indicates that there is a time delay between short 
burst occurrence and the main epoch of star formation 
activity in the hosts, as determined by the progenitor's 
lifetime. 

In this paper, we present a method to deduce the life- 
time of short GRB progenitors from the burst rates in 
host galaxies of different types that have distinct star 
formation histories. Limited by the accuracy of how well 
we can separate their contribution to the star formation 
history, we show, with this constraint in mind, how a 
large sample host galaxies could be used to determine 
the lifetime distribution of short GRB progenitors. Here 
we chose to separate the formation history of stars re- 
siding in early- and late- type galaxies, but, as shown in 
§5, our study can be easily generalized to any set of sub- 
populations of galaxies possessing distinct star formation 
histories. The layout is as follows. In § 2, we review 
the local stellar budgets in early- and late- type galaxies. 
The total star formation history is subsequently decom- 
posed in § 3 into a sum of the stars that assembled in 
today's early- and late- type galaxies. With the decom- 
position in place, in § 4, we show how the frequency of 
short GRBs in a well-defined sample of host galaxies of 
different types can be used to constrain the lifetime of 
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Fig. 1. — The local stellar mass density from galaxies of different 
types. The histogram with errorbars is for galaxies of all types 
while the dashed curve is a Schechter function fit to the stellar mass 
function of the entire galaxy sample. Red and blue histograms are 
for early- and late-type galaxies, respectively. The early- and late- 
type galaxies are divided according to their color (dotted) or light 
concentration (solid) . This plot is bas ed on the galaxy stellar mass 
functions derived bv lBell et all 112001) . 

the progenitors. Finally, in § 5 we summarize our re- 
sults and outline future prospects. Throughout this pa- 
per, we assume a spatially-flat ACDM cosmology with 
Q m = 1 — Ha = 0.3 and the Hubble constant h = 0.7 in 
units of lOOkms^Mpc -1 . 

2. THE LOCAL STELLAR MASS INVENTORY 

The focus of this paper is the relatively local (z ~ 0) 
host galaxy population, but our study can be easily gen- 
eralized to any redshift. The main idea is to use the 
difference in the star formation history of galaxy sub- 
populations to probe the lifetime distribution of short 
GRB progenitors. As an example, here we use the event 
rate of burst residing in local early- and late-type galax- 
ies. On average, early-type galaxies have their stars 
formed earlier than late-type galaxies, and this differ- 
ence, together with the delay time for short GRB out- 
burst, leads to different burst rates in these two types 
of galaxies. As shown in § 3, our method of inferring 
the star formation history of a given population relies 
on tracing the assembly's history of the stellar mass in 
galaxies at z = 0. For this reason, in this section, we 
briefly review the local stellar mass budgets in galaxies 
of different types. 

With a universally applicable stellar initial mass func- 
tion (IMF), the stellar mass function (MF) of galaxies 
can be estimated from well-defined samples of galax- 
ies. Galaxy MFs have been meas ured based on a 
few l a rge galaxy redsh ift surveys (e.g.. [Kauffman n et al.l 
120031: iBell et all koiKSli. For the calculations presented 
here, we adopt the galaxy MF measured by IBell et afl 
(2003) using a large sample of galaxies from the Two Mi- 
cron All Sky Survey (2MASS; ISkrutskie et"aU[2006l) an d 
the Sloan Digital Sky Survey fSDSS: lYork et all 120001 1. 
The Bell et al. catalog provides well-defined samples of 
galaxies of different types as defined by either the light 
concentration or the color of galaxies. 



In Figure [l] we show the stellar MF, <j>, multiplied by 
the stellar mass, M star , for galaxies of different types. 
The quantity p s %ax{M s%ax ) = M siaI ^(M sta . r ) is tne stellar 
mass density contributed by galaxies with stellar mass 
Mstar- The peak of p st ar for galaxies of all types comes 
from M* galaxies, where M* ~ lO n M0 is the charac- 
teristic stellar mass derived from the Schechter function 
fit. The stellar mass density distribution for late- type 
galaxies peaks at masses slightly less than M* and has a 
tail extending to low stellar masses. Early-type galaxies 
have a narrow stellar mass density distribution around 
M* and dominate the local stellar mass budget above a 
few times 1O 1O M0. Although early-type galaxies are far 
less in total number than late- type galaxies, they are on 
average more massive and thus make a larger contribu- 
tion to the local stellar mass density. The total stellar 
mass from early-type galaxies is about 1.3 (2.3) times 
that from late-type galaxies, if the classification of galaxy 
types is based on light concentration (color). 

If early- and late-type galaxies had similar star forma- 
tion histories, then the ratio of the short GRB rates in 
galaxies of different types would simply be given by the 
ratio of the stellar mass density in these galaxies. How- 
ever, the two types of galaxies have distinct star forma- 
tion histories - on average, stellar populations in early- 
type galaxies are older than those in late- type galaxies. 
For this reason, the delay time between formation and 
the short GRB outburst plays an important role in de- 
termining the burst rates in these two types of galaxies. 

3. THE STAR FORMATION HISTORY OF GALAXIES 

The cosmic star formation history (SFH), which in- 
cludes stars forming in both early- and late-type galax- 
ies, can be probed by different techniques, including the 
rest-frame UV continuum emitted by young stars, the in- 
frared or sub-millimeter reprocessed radia tion, and line 
emission from star-forming regions (e.g.. iMadau et 
1 19961: iGiavalisco et al.ll2004t iPerez-Gonzalez et al.ll2005t 
Hopkins et all 120001 : see lHopkinsll2004l and lFardal et all 
120061 for a summary and a compilation of SFHs from 
different techniques). 

Analytic fits to the SFH based on compiled star forma- 
tion rate (SFR) at va rious redshifts have been derived by 
many authors (e.g., E ppkins 2004: Hop kins fc Beacoml 
1200a ICole et all l200lTiFardal et al.ll2006fl . Here, as an 
illustrati on, we use the ana lytic fit to the overall SFH 
given bv lFardal et~aTl (|2006h . 



SFR aU (z) 



PlP2P3Po(l + z) 1 



[l+Pl(l + z)P2]P3 



(1) 



where p = 9.0 x 10 8 M Q Mpc _d , p l = 0.075, p 2 = 3.7, 
P3 = 0.84, and H(z) is the Hubble constant at redshift z. 
The above fit assumes a diet Salpeter stellar initial mass 
function 4 (IMF), To the first order, altering the IMF 
leads to an overall vertical shift of the SFR. For example, 
adopting a Salpeter IMF, the SFR would be about 40% 
higher. What is more, even with a fixed IMF, the derived 
SFR as a function of redshift shows considerable scatter. 
As a result, the SFR in equ ation (U) may be s caled up 
and down by about 0.2 dex (|Fardalet alj2006"t ). 



4 A diet Salpeter IMF is similar to the typical Salpeter IMF at 
the high mass end but differs at the low end by having less lower 
mass stars so that t he total stellar mass c ontent is 0.7 times that 
of the Salpeter IMF (IBell fc de Jondl200lh . 
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In addition to the total SFH, we are also interested in 
knowing the SFH for stars residing in galaxies of different 
types. More specifically, given the z ~ early- and late- 
type galaxy populations, our aim is to understand how 
the stars in these two different types of galaxies were 
pieced together. We, of course, only need to know the 
SFH for galaxies of one type, as the SFH for the other 
type can be easily derived by subtracting the known con- 
tribution from the total SFH in equation |T]). 

Ideally, to determine the SFH of a given population 
of z ~ galaxies, one would like to empirically derive 
their assembly history by identifying their progenitor sys- 
tems at different redshifts. Alternatively, we can study 
the stellar population contents of z ~ galaxies to infer 
their SFH. The former is not trivial, while considerable 
efforts have been mad e to the latter through spectral syn- 
thesi s modeling (e.g.. IHeavens et alj|2004t iPanter et al.l 
I2006t ICid Fernandes et all 120051 I2006D . which decom- 
poses the observed spectra of galaxies into stellar popula- 
tions of a range of ages and metallicities. iHeavens et alJ 
(|2004f ) first used this technique to infer the global SFH 
and that as a func tion of stellar mass. More recently, 
IPanter et al.l (|2006f ) improved it to yield a better agree- 
ment with other empirical determinations of the global 
SFH, w hich peaks at systema tically higher redshifts than 
that in IHeavens et al.l (|2004l ). Obviously, this method 
would improve with a better understanding of the sys- 
tematic effects. 

A more readily available method is to use galaxy for- 
mation models to separate the SFH as a function of 
galaxy type. The accuracy of galaxy formation models 
crucially depends on our understanding of the gas physics 
and star formation as well as feedback processes. Cur- 
rent galaxy formation models, although far from provid- 
ing a perfect description, can successfully explain many 
aspects of the observed properties of galaxies. Hope- 
fully, future improvements in galaxy formation modeling 
would help provide a better description of the SFH as a 
function of galaxy type. 

With this in mind, we proceed separating the global 
SFH using a particular galaxy formation model (with 
model uncertainties roughly accounted for). It should 
be noted here that our method is applicable to any two 
or more galaxy populations provided they have distinct 
SFHs. This is clearly illustrated in § 5, in which, for 
comparison, an observationally motivated decomposition 
b ased on galaxy masses is used instead. 

iDe Lucia et al.l ([2006) studied the formation history of 
elliptical galaxies using a galaxy formation model based 
on the Mille nnium Simulation of the concordance ACDM 
cosmology (jSpringel et alJ 12005). They calculated the 
average SFH of z ~ elliptical galaxies of various stellar 
masses. If early-type galaxies arc identified here as ellip- 
ticals, we can then compute their average star formation 
rate (M0yr _1 Mpc~ 3 ) as a function of redshift using 

SFR carly (z) = / rfM atar dF ^ Stan Z) M star <MM star ), 

J dMstar 

(2) 

where 0B(M st ar) is the z = ste ll ar M F of 
early-type galaxies derived by Be ll et all ([2003!), and 
dF(M Btal , z)/dM star (in units of A^yr" 1 !^) is the av- 
erage SFR per stellar mass for a z = elliptical galaxy 
of stellar mass M sta r (obtained from spline interpolations 



between the curves presented in lDe Lucia et alj |2006). 

The application of equation ([2|) combines the empiri- 
cally inferred galaxy stellar mass function and the the- 
oretically predicted SFR. It requires the galaxy samples 
derived by both observation and theory to be consis- 
tent. Although in principle this is not a problem, since 
any sophisticated galaxy formation model would give de- 
tail predictions for a well defined galaxy population, the 
available data and models used here are based on slightly 
different definition rules, which introduces uncertainties 
in the determinatio n of the SFH of early-type galaxies. 
In lBell et alJ (|2003h , galaxies are classified through either 
color or light concentration, which is shown to be only 
a cr ude discriminant bet ween early- and late-type galax- 
ies (IStrateya et alj|2001f>. an d, as a result, the early- type 



MM) include galaxies earlier than 
20061 ) split galaxies in their forma- 



gal axies in Bell et al 
Sa. IDe Lucia et al.l 
tion model according to the predicted B-band bulge-to- 
disk light ratio and classify them as ellipticals provided 
the bulge light is at least 70% of the total light. For this 
reason, the SFH of the observed local early-type galaxies 
as derived from the stellar-mass dependent SFH of the 
modeled ellipticals would tend to shift toward slightly 
earlier formation times. In addition, in the galaxy for- 
mation model used here, stars are produced from cooled 
gas and a stellar IMF is adopted only for the purpose of 
providing spectroscopic predictions. Because of this, a 
comparison between this IMF and the one used to derive 
the observed galaxy stellar mass function is not mean- 
ingful. Instead, the key questions are whether the IMF 
assumed to derive the observed galaxy mass function is 
close to the "true" one and whether the star formation 
efficiency in the galaxy formation model is accurate. In- 
accuracy in either one would lead to a relative amplitude 
change in the SFHs. 

The left panel of Figure [2] shows an example of the SFH 
decomposition. Although the primary aim of the SFH 
decomposition selected in this study is to illustrate the 
applicability of our method, we pay particular attention 
to the uncertainties associated with it. Several variants 
of the SFH decomposition are constructed to test the ro- 
bustness of the results. Each version has a different way 
to adjust the derived SFR in early-type galaxies with re- 
spect to the assumed overall SFH. The amplitude and/or 
position of the derived rate is adjusted using one of the 
following methods: (i) the derived SFR is selected to be 
equal to the overall SFH whenever the predicted rate lies 
above it (red solid curve in Fig. [2]) , (ii) the overall am- 
plitude of the derived SFR is systematically lowered so 
that it lies just below the overall SFH (red dotted curve 
in Fig. [2]), and (iii) the derived SFR is arbitrarily shifted 
to low redshifts so that it both lies below the overall 
SFH and peaks at a redshift that is not lower than the 
overall SFH (red dashed curve in Fig. [2]). In addition, 
the overall SFH is allowed to change by 0.2 dex in order 
to account f or the uncertaintie s associated with its de- 
termination ([Fardal et a l. 2006). The above variants are 
applied every-time the overall SFH is altered. 

Finally, we adopt the tw o definitions of early-type 
galaxies in iBell et ah! (|2003f ) based on color and concen- 
tration, respectively. As a whole, the above procedure 
takes into account both the inaccuracy of the galaxy for- 
mation model predictions and the uncertainties in the 
global SFH as inferred from observations, as well as the 
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Fig. 2. — Left panel: Star formation histories for z ~ galaxies of different types. The total star formation history (black) is decomposed 
into two parts — formation histories of stars that are now in the z ~ early- (red) and late-type galaxies (blue). Solid, dotted and dashed 
red and blue curves result from different treatments of the SFH of early-type galaxies (see text), which give a clear representation of the 
uncertainties associated with the assumed SFH decomposition. The curves shown here are for early- and late-type galaxies as defined by 
the concentration of their light profile (corresponding to the solid histograms in FigUJ. Right panel: The ratio of short GRB rates in early- 
and late-type galaxies at z ~ as a function of the index n of the progenitor lifetime distribution P(t) oc t™. The thick solid, dotted 
and dashed curves correspond to the different SFH splits shown in the left panel. The thin curves, on the other hand, are for the same 
SFH decomposition but for galaxy types classified based on color instead of light concentration. The two dot-dashed curves give a fair 
representation of the uncertainties in the determination of the burst ratio, which arise primarily from the different methods and adjustments 
used here to decompose the SFH (see text). 



ambiguity of the earl-type definition criteria. These un- 
certainties are propagated self consistently when estimat- 
ing the lifetime distribution of short GRB progenitors 
(§HJ. The region enclosed by the dashed-dotted curves 
in the right panel of Fig. [2] shows the combined effects on 
the derived lifetime distributions of short GRB progen- 
itors when adopting the three variants discussed above, 
allowing a 0.2 dex change in the overall SFH, and con- 
sidering the two definitions of early-type galaxies. 

As a representative example, in the left panel of Fig- 
ureEl we plot the SFH in equation (JTJ) and the SFH of the 
z = early-type galaxies defined through color, as black 
and red curves, respectively. The different line types for 
the SFH of early-type galaxies are for the various ad- 
justments listed above. In general, the derived average 
SFHs of the z = early-type galaxies peak at a look-back 
time of about ~10-12 Gyr and subsequently decline to- 
ward low redshifts. The shape of this early-type SFH is 
qu ite consistent wi t h tha t of the fossil bulge model fit 
in INagamine et al.1 (|2006l ) and that of SF H of galaxies 
with s tellar mass above 1O 11 M inferred in lPanter et al.l 

The average SFH for the z = late-type galaxies, 
computed by subtracting © from (Q]), is shown as blue 
curves in the left panel of Figure [5] The average SFR for 
the z — late-type galaxies has a steep rising starting 
at lookback times of ~11 Gyr, subsequently peaking at 
~8-9 Gyr, and slowly declining thereafter. The SFR be- 
comes increasingly dominated by late-type galaxies from 
lookback times of ~9 Gyr to the present epoch. 

We note here that the overall SFH in equation (pTJ) has 
a slightly higher amplitude than the predicted early-type 
SFR at lookback times greater than 12 Gyr, which gives 
rise to a modest second peak in the late-type SFR at high 
redshifts. This may be interpreted as a contribution from 



local late-type galaxies to the SFR at these early epochs, 
which could be related to the formation of the bulges 
of these galaxies. It may also be an artifact associated 
with the inaccuracy of the model, which we consider to 
be more likely. We find that even if taken at face value, 
this early contribution to the formation of stars in local 
late-type galaxies leads only to a less than 10% increase 
in the inferred short GRB burst ratio in the two types 
of galaxies (§ |4|). For this reason, we leave the SFR for 
late-type galaxies at these epochs as those given by the 
decomposition and do not perform any correction. 

4. CONSTRAINTS ON THE LIFETIME OF SHORT GRB 
PROGENITORS 

The short GRB rate is a convolution of the SFR and 
the distributio n of time d elays between formation and 
outburst fe.g.. lPiranlll992D 



Ri{z) = C 



t(z) 



dr SFRj(t - t)P(t), 



(3) 



where the subscript i denotes the type of galaxies in 
consideration, P{t) is the probability distribution of the 
time delay r, and C is a normalization constant. In prin- 
ciple, details of the star formation process may be related 
to the assembly history of galaxies. Therefore, the dis- 
tribution P(t) and the normalization constant C could 
be different for different types of galaxies. We make the 
simple assumption here that this dependence is weak so 
that the same distribution and normalization is used for 
all types of galaxies. 

The distribution P(t) of the time delay for short GRB 
is not yet well understood. It can be constrained us- 
ing, for e xample, the burst rate as a function of red- 
shift (e.g., iNakar et al.ll200llGuetta fc Piranll2005l:lAndol 
120041 ). The luminosity function of short GRBs and its 
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redshift evolution used in this method are not well deter- 
mined at present. To avoid such complications, here we 
propose to use the local rates of short GRBs in different 
types of galaxies as an alternative and complementary 
method to constrain the probability distribution P(t). 
As an illustration, we adopt a simple parameterization, 
P(r) oc t", and calculate the ratio of rates of short GRBs 
at z ~ in early- and late-type galaxies as a function of 
n. The right panel of Figure [2] shows the results of such 
exercise. The thick solid, dotted, and dashed curves are 
for three different treatments of the early-type SFR (sec 
§ E]) , corresponding to the three different curves in the 
left panel. Here galaxy types are defined by their light 
concentration. The three thin curves in the right panel 
have been calculated using a similar treatment but are 
for early-type galaxies as defined by their color. The 
two dot-dashed curves give the maximum range of un- 
certainties in the resultant burst rate ratio, as they give 
a fair representation of all the uncertainties involved in 
separating the SFH as discussed in § 3. 

For larger n, the distribution P(r) is weighted more to- 
ward longer time delays, so that early-type z ~ galax- 
ies, which on average form their stars earlier than late- 
type z ~ galaxies, are more likely to host short GRBs. 
For smaller (negative) n, short GRB progenitor systems 
with shorter time delays would dominate and one would 
be more likely to find them in late-type galaxies. For 
n < —3/2, short GRB progenitors would be those with 
very short time delays so that the short GRB rate ratio 
is basically determined by the ratio of the z ~ SFRs in 
early- and late-type galaxies. However, the amplitude of 
the plateau at n < —3/2 in the right panel of Figure [2] 
should be taken cautiously, because it depends on how 
on e interpolates the SFR in early-type galaxies derived 
bv lDe Lucia et alj (|2006l ) to z ~ 0. For this reason, the 
method proposed here becomes less useful for time delays 
that are short when comparing to the mean age of stars 
residing in early- type galaxies. Clearly, as we argue in 
what follows, current observational constraints for short 
GRBs are such that uncertainties in the determination of 
the amplitude of the plateau at n < —3/2 are irrelevant. 

Our calculation is based on z ~ galaxy populations. 
Of the nine well localized short bursts, only three hosts, 
two elliptical galaxies and one star-forming galaxy, re- 
side at z < 0.3. Therefore, the face value of the burst 
ratio in local early- and late-type galaxies is ~ 2 provided 
there are no selection effects that make short GRBs more 
likely to be detected in a galaxy of a given type (e.g., 
because of the low density of the interstellar medium 
one may expect the GRB afterglow in e lliptical galaxies 
to be faint, e.g.. iBelczvnski et al.ll200"6f ). The error-bar 
on this observed ratio is not readily assigned because 
of the small number statistics. However, even with a 
conservative estimate of the ratio as ~ 2 ± 1.5, the anal- 
ysis presented here alludes to a probability distribution 
with n > — 1, thus favoring long delay times. This is 
in contrast to what is found in Type la supernovae, for 
which the relatively low fr e quency in early-type g alaxies 
(e.g., Ivan den Berehlll990t iMannucci etafl I2005D yields 
n < -1. 

If the distribution P(t) is assumed to have a log- 
normal form, there are a series of combinations of the 
mean and the width that can explain the observed burst 
rate ratio in early- and late- type galaxies at z ~ 0. We 
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Fig. 3. — Illustration of different functional forms of P(r) that 
give rise to the same observed short burst rate ratio in early- and 
late-type galaxies. If P(t) is assumed to have a log-normal form, 
there are a series of combinations of the mean and the width that 
can explain a given observed burst rate ratio, here assumed to be 
either 1 (bottom panel) or 2 (top panel). A narrow and a wide 
log-normal distributions are shown. 

illustrate the underlying powerlaw and two cases of log- 
normal P(t) distributions for burst ratio of one and two, 
respectively, in Figure [3l using the SFH split correspond- 
ing to the top curve in the right panel of Figure [U For a 
narrow log-normal distribution, it just peaks at tq when 
the ratio of the SFRs in two types of galaxies at a look- 
back time To equals the burst ratio. A wide log-normal 
distribution and the corresponding powerlaw distribution 
have a general similarity with each other, but they di- 
verge from each other at small time delays if the power- 
law index is negative, which indicates that the method 
has a poor constraint in P(t) at the low r end. For a 
burst ratio of unity, we find n > —0.45. If we take 13.7 
Gyr as an upper limit for the delay time, which allows 
short GRB progenitors to outburst within a Hubble time, 
a distribution P(r) with n = —0.45 implies that about 
half of progenitors of short GRBs should have lifetimes 
longer than 3.9Gyr. If at lower r, P(r) does not rise as 
steeply as a powerlaw, this number becomes higher — for 
example, with a cutoff in P(r) at r <2Gyr, the number 
becomes 6.7Gyr. For a wide (narrow) log-normal distri- 
bution, a burst ratio of order unity implies that about 
half of progenitors of short GRBs should have lifetimes 
longer than 7.0 (8.5) Gyr. 

Obviously, the above calculation is only sketchy and 
should be taken with a grain of salt at present, al- 
though, given the uncertainties assigned to the modeling, 
the lower bound of n seems rather robust. Constraints 
should improve as more host galaxies of short GRBs are 
detected and the SFH separation method is improved. 
In order to reach a 10% accuracy in the the burst ratio, 
observations of some 100 local short GRB host galaxies 
are needed. If the split of the SFHs could be well deter- 
mined, an associated uncertainty An ~ 0.3 in the time 
delay distribution could be easily achieved. If the real- 
istic SFH separation lead to a 20% uncertainty in the 
predicted burst rate ratio, the constraints on n would be 
degraded to An ~ 0.9. On the other hand, as suggested 
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by our analysis, it is possible to set a very reliable lower 
bound on n. In brief, for a more precise constraint to be 
derived on the time delay of short GRB progenitors using 
the method proposed here, one needs a larger sample of 
local short GRB host galaxies to be observed, the split of 
the SFH to be better determined, and the classification 
scheme of GRB host galaxies to be consistent with the 
assumed SFH split. 

5. DISCUSSION 

In the local universe, based on the stellar mass func- 
tions of galaxies, about 55-70% of the stellar mass is 
in early-type galaxies, and, from the SFH we inferred, 
the corresponding stars mainly formed about 9 Gyr ago. 
Two of the three host galaxies of local short GRBs 
(z < 0.3) are associated with old and massive galax- 
ies with little current or recent star formation, which 
makes it unlikely that short bursts are associated with 
massive stars. Presently available data suggests, but not 
yet prove, a long time delay between the formation of 
the progenitor system and the short GRB outburst - 
for progenitors that can outburst within a Hubble time, 
about half of them have lifetime longer than ~7 Gyr. It 
is fair to conclude that, based on the current host galaxy 
sample, the progenitors of short GRBs app ears to be 
longer lived than those of Type la supernovae. I Fox et all 
(2005) also reach similar conclusion by arguing that Type 
la supernovae occur more frequently in late-type, star- 
forming galaxies. Simply based on a comparison of Hub- 
ble types bet ween short GRB s and Type la supernovae 
host galaxies, iGal-Yam et al.l (|2005h argue that the de- 
lay time of short GRBs should be several Gyr even if the 
Type la supernovae delay time is as short as ~ 1 Gyr. 

The lifetime of the progenitor systems is estimated here 
by using the SFH of elliptical galaxies from a galaxy for- 
mation model. This allow us to separate the early- and 
late-type galaxy contributions to the overall cosmic SFH. 
It would, however, be more self-consistent to infer SFHs 
of different types of galaxies by modeling the observed 
spectra with stellar population synthesis models. In ei- 
ther case, the uncertainty in the derived SFHs should 
be folded into the errors derived by this method for the 
distribution of time delays of short GRB progenitors. 

In our calculation, different definitions of early-type 
galaxies (by color or by light concentration) introduce 
uncertainties in the resultant SFH. This, in principle, 
would not be a problem since we can choose to use the 
same definition for identifying the short GRB host galaxy 
type. The method proposed here is not limited to a sepa- 
ration between early and late galaxies. As long as galax- 
ies are divided into two (or more) sub-populations that 
have distinct star formation histories, they can be used 
to constrain the lifetime distribution of GRB (or Type 
la supernova) progenitors, foll owing the same rea soning 
presented here. For example, IShin fc Bergerl (|2006[ ) re- 
cently applied our approach to field and cluster elliptical 
galaxies. Another example would be to divide galax- 
ies according to their stellar mass. iPanter et al.l (2006) 
present SFH of local galaxies as a function of galaxy stel- 
lar mass by modeling their spectra. Low mass galaxies 
on average form their stars later than high mass galaxies, 
which can be used to constrain P(t) if the stellar mass of 
GRB hosts can be obtained. As an illustration, in Fig- 
ure 21 we show the separation of the SFR of galaxies with 



low and high s tellar masses (belo w or above 10 11 M Q ) us- 
ing results bv lPanter et al.l (l2006h. T he shape of the over- 
all SFH used in IPanter etal . (2006) is slightly different 
from that in Figure [2j which may reflect some systemat- 
ics in the derivation of the SFH based on the fossil record. 
If the systematics can be well controlled, this approach 
can become even more powerful than the one used here 
by inferring SFHs of individual GRB host galaxies. A 
maximum likelihood method can then be used to con- 
strain P(t) based solely on the SFH of host galaxies, 
which puts our proposed method to an extreme — divid- 
ing the galaxy populations into individual galaxies. Fur- 
thermore, instead of assuming a functional form of P{t), 
such an application to a large number of individual host 
galaxies may allow constraints on a non-parametric form 
of the distribution [i.e., constraining P(r) in different r 
bins] . 

In this paper, we limit our study to z ~ galaxies, 
but the method can be easily generalized to galaxy pop- 
ulations at any redshift provided that one can accurately 
infer their SFHs (e.g., through spectral synthesis). The 
observational task can be minimized by focusing on the 
observations of host galaxies of individual short GRBs. 
However, for applications at high redshift to be useful, 
the luminosity function of short GRBs needs to be un- 
derstood. 

Throughout this paper we have assumed the same life- 
time distribution of short GRB progenitors in both early- 
and late-type galaxies. However, star formation pro- 
cesses in these two types of galaxies may not be identical. 
For example, elliptical galaxi es can form by the mergin g 
of two gas-rich galaxies (e.g., iMihos &; Hernquistl l'1994). 
Many globular c lusters can form in the merging process 
(jSchwe izcr 2003), which could enhance, for example, the 
fraction of binary pro genitors and also cha nge the life- 
time distribution (e.g.. iGrindlav et al.ll2006l ). The mag- 
nitude of this kind of effect on the GRB progenitors is 
a formidable challenge to theorists and to computational 
techniques. It is, also, a formidable challenge for ob- 
servers, in their quest for detecting minute details in ex- 
tremely faint and distant sources. 

At least two short GRB host galaxies are found in clus- 
ter environments. There may exist a sele ction bias of 
detec ting short GRBs in a dense medium (jBloom et al.1 
2006). To study the association of short GRBs with 
clusters, it would be useful to separate the stellar mass 
function into that for field galaxies and that for clus- 
ter galaxies in addition to early- and late- type galaxies. 
More promising for the immediate future, the prepon- 
derance of cluster environments can be investigated ob- 
servationally. Important information may be gained by 
studying the local stellar mass inventory shown in Fig- 
ure [TJ Approximately 50% of the stellar mass contents in 
early-type galaxies are in galaxies with M sta r > 10 11 M Q 
that typically reside in clusters. Since it is likely that 
in clusters galaxies shut off their star formation process 
early on, a long progenitor lifetime further increases the 
tendency for short GRBs to happen in cluster galaxies. 
It is fair to conclude that the observed preponderance of 
cluster environments for short GRBs is consistent with 
an old stellar population that preferentially resides in 
early- type galaxies. 

Detailed observations of the astrophysics of individual 
GRB host galaxies may be essential before stringent con- 
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Fig. 4. — Similar t o Fig, [2] but the S FH is decomposed into those for galaxies with different stellar masses. This decomposition is 
based on the results in Pantcr et al. (2006), where the SFHs of local, individual galaxies are inferred from modeling their spectra and are 
then subsequent ly grouped accordin g to the present stellar mass content. In the left panel, the solid curve shows a spline fit to the overall 
SFH inferred bv lPanter et al.l 1)20061 ) . and the dashed and dotted curves give the spline fits to the SFRs of galaxies with stellar mass below 
and above 10 11 Mq, respectively. For lookback time larger than ~12.5 Gyr, a cutoff profile is introduced. In the right panel, the ratio of 
short GRB rates in galaxies with stellar mass above and below 10 11 Mq is plotted as a function of the index n of the progenitor lifetime 
distribution P(t) oc t™. The purpose of the figure is to illustrate that the method presented in this paper can be generalized to any 
sub-populations of galaxies that have distinct SFHs (see text). 



straints on the lifetime of short GRB progenitors can be 
placed. If confirmed with further host observations, this 
tendency of short GRB progenitors to be relatively old 
can help differentiate between various ways of forming a 
short GRB. 
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